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ABSTRACT 


A theory is developed which describes the local wave 
effects near a wave front advancing into calm water and is 
uniformly valid near the front. The waves generated by an 
oscillating vertical boundary (a wavemaker) and propagating 
out into an infinitely long channel are investigated both 
theoretically and experimentally. Comparison of the theo- 
retical and experimental results indicates that the theory 
does predict the wave phenomena about the front to a good 
degree of accuracy. 


Conclusions are drawn regarding the transient wave 
response to a wavemaker and compared to those of John W. 
Miles regarding the transient gravity wave response to an 
oscillating pressure. The transient features of the wave 
front considered are the width of the transient zone at a 
given time since the disturbance to the free surface of the 
fluid and the length of the transient period at a given 
distance from the disturbance. It is shown that the con- 
clusions regarding transient gravity wave phenomena apply 
without qualification to the waves generated by both oscil- 
lating pressure sources and moving surfaces. 
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NOMENCLATURE 


6- Dirac @deltaviunction 


n- Wave height 


g - Acceleration due to gravity 


G - Green's function 


r- Gamma function 


h 


H 


depth of fluid 

Heaviside step function 

Imaginary part of a complex variable 
Real part of a complex variable 
Displacement of wavemaker paddle 


Amplitude of displacement of wavemaker 
paddle at free surface 


The ratio of distance to time (x/t) 
Same as V above. 


Frequency of oscillation of wavemaker paddle 


j~ 
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INTRODUCTION 


Tne object of this study is to develop an expression 
Winicn describes the local wave effects near a wave front 
aavancing into calmowater thateis uniformly valid aboutritne 
front. Limitime formswef this expression aS — > © or x > © 
were aeveloped nistorically by Lamb~ and by Cauchy and 
Poisson respectively where t represents elapsed time after 
an initial disturbance to tne free surface of the fluid and 
x the distance from the disturbance. Lamb's result de- 
scribes tne waves resulting from a steadily oscillating 
disturbance where the elapsed time since the disturbance 
was initiated is assumed to be so large that all transient 
features have disappeared ana only tne so-called steady 
state phenomena remains. The Cauchy-Polsson expression 
describes the initial propagation of waves out from a dis- 
turbance but does not account for the time dependent aspects 
of tne disturbance. Tne exact formulation of these expres- 
Sions depends on the nature of the disturbance, e.g., 
pressure variations or moving surface - and are widely 


Puodlable in the literature. 


1 


Lamb, H., "On Deep Water Waves", Proc. Lond. Math. Soe., 
5-2, 371,463). 
“caucity, A., Pewmsson, S. D. "Memoire sur al theorie des 


ondes", Mem. de Ll'Acad. Roy. des Setences, i. (1816). 
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Neither of the above theories attempt to describe the 
transient wave phenomena spanning their respective limits. 
This study develops a theory that predicts the waves in 
this transient zone and which upon passage to the appro- 
priate limit yields the nistorical result. The waves gen- 
erated by an oscillating vertical boundary (a wavemaker) 
and propagating out into an infinitely long channel are 
investigated both theoretically and experimentally. The 
channel is assumed to be infinitely wide so that the 


problem is two-dimensional. 


During the problem formulation phase of this analysis, 
tne attention of the author was directed to a paper by 
Miles”, which considered the transient gravity wave response 
to an oscillating pressure. Although the mechanism by which 
transient waves are generated is different in Miles' 
approach, the qualitative features of the resulting waves 
snoulu be similar to tnose generated by a wavemaker. In nis 
paper Miles developed an expression predicting the transient 
waves thac was uniformly valid about the wavefront. His 
solution method, though apparently consistent, was not 
immediately satisfying as it involved finding the asymptotic 
expansion of a double integral by applying the method of 


stationary phase twice in succession. It was felt that a 


“Miles, J. W., "Transient Gravity Wave Response to an Oscil- 
faerie resscUromemin 8 77 Vech., 13, (1962), pp. T45sloue 








more direct means of obtaining a uniformly valid expres- 
sion was desirable - both to verify Miles' result and to 


obtain clearer insight into the transient problem. 


Tnis study was undertaken primarily in the interest of 
academic completeness; viz., to bridge the gap between the 
solutions of Lamb and Cauchy-Poisson. tHWowever, several 
engineering applications can be foreseen. The utilization 
of wavemakers in towing tanks with corresponding applica- 
tions such as ship model and ocean structure testing is of 
direct interest to the naval architect or ocean engineer. 
Transient features of generated waves may allow some 
aspects of model testing to be explored more rigorously and 
completely. Further application is found in the area of 
ocean wave forecasting. The wave generating models used 
presently in predicting the wave phenomena associated with 
a storm at sea are not completely satisfactory and a new 
model, based partially or entirely on the expression 


derived in this study could be formulated. 





THEORETICAL ANALYSIS 


Consider the waves generated by a vertical plate in 
oscillatory motion given by (1) in a fluid bounded by a 
free surface, tne vertical plate, and of infinite extent 
Otherwise. Tinis so-called wavemaker problem is solved in 
Appendix 2B and the resulting expression for tne waveheignt 


as a function of distance and time is given by (2). 


(1) S. <ic e a ie oreni(ee mECen 


wnere S = amplitude of displacement at the 
free surface 


H(t) is the Heaviside step function 


(2) re _ wS ~ cos (ux) sin(wt) - sin(/Ygu t) : 
o U ti 7g w - vYgu 


i Sin(wt) + sSin(vgu t) an 
w + Ygu 
Introducing the change in variable u = o7/g 
ee a a } oon 
ey Gono = 2ws SCOS| (Gn 7G) <4 sin(wt) sin(ot) - 
T o* + w* w- 0 
O 
i‘ Sin(wt) + sin(ot) sais 
Weed 0 
ww 2 2 2 ® 
i(@xtut) i (2x-ot) i (Lx+wt) i Cx-wt) 
_1 {2wS| o2e mee =o ae 
mj) 7 


(o7+w*) (otw) (o-w) 





Using nondimensional variables x = w*x/g , t =wt, ~ 
me O/w and nh = x , (3) reduces to 
(4) 
OR ey) ae 


iG 





« 2 ‘ e 2 « 2 . 2 
+ = —_ 
: os xt+ot) ae ito -xX—oe) 3 1 UC t } 5 L(6 "x= 
m 


‘ (o*+1) (o+1) (0-1) 


Equation (4) for n(x,t) may be considered as the sum 


of four integrals. 











(5) i (xa I (I, - I, = I, + ae 
CO : 2 
where I, = | £ Co) arid RTOu) 
O 
I, _ | een iO - x Ot) 5 
O 
2 . 2 
I, = | sh pen he te) 
O 
co é 2 = 
7) Z | = {91 ai ld x t) ag 
O 
; cE . 
ancgeweet (co) = 
go t=] 


The asymptotic behavior of these integrals for 
large x or t 1s readily evaluated by treating 6 as -ascom= 
plex variable and employing tne residue theorem and 
Cauchy's integral theorem of complex variable theory. 
Their behavior for small x or t is given in Appendix H. 
In the coinplex o plane, the real integrals in (5) 


are transformed to complex integrals which are to be 





evaluated along the contour consisting of the ray from the 


Origin to infinity along the positive real axis. 


It is easily shown (see Appendix C) that the asymptotic 
behavior of integrals Ij, 1, and I, is given by the follow- 
ing expressions for large t or x: 


4 a ai (xtt) + 0(t73) 


ee El el (Xtt) 4 9 (4-2) 


oe ei (x-t) + 0(t72) 


Tne asymptotic behavior of (5) then reduces to the 


following: 


(6) n(x,t) v IL at oe > is \ 


The asymptotic behavior of I, is not as easily as- 


2 
1’ I, and I,- In the evaluation 


of these integrals, the best choice of an equivalent con- 


certained as was that of I 


tour to facilitate the determination of their asymptotic 
behavior was obvious. If a similar contour is chosen for 
tne evaluation of I., however, tne integral diverges along 
portions of the contour. Consider then the exponential 
term of I. in expanded form where the substitution 

ag 
e 


se- © has been maae. It takes the form 





(7) ei (o*x-ot) = aitir® (x/t) cos (26) -r sin(6@)} x 


x qe tix® (x/t)sin(26)—xr sin(6) } 


A good choice for an equivalent contour would be one on 
which this term becomes small as x and t become large and 
vanishes when the contour is extended to infinity. For 
this term to become exponentially small for large values of 


t (or x), the following restriction on r and @ applies: 


(8) Y-(x7e) Simeee ye = Sinte) > 0 


“A 


Let x/t = V: then for 6 negative r 1/2V cos (86) 


for 6 positive r > 1/2V cos(6) 


for 6 = 0 iG Wy 


Thus, any contour with the general features of Co 
shown below 1S a satisfactory equivalent contour for the 
Pac from 0 to + © along the real axis. The point 


o = 1/2V + Oi where the contour C. crosses the real axis is 


Z 
a stationary point of the integrand and must correspond 
to a saddle point as the integrand is a function of the 


complex variable o and therefore satisfies Laplace's 


equation. 





o-plane 





Here Singularities of tne integrand of I, have been 
indicated with x marks and in all cases are simple poles. 
The contour along the real axis is denoted by C, and has 
been indented at o = 1. Thus, the above integral is 
defined in the Cauchy principle value sense. If the 
integral exists in the conventional sense, then its value 


will necessarily agree with the Cauchy principle value so 


there is no loss in generality. 


The existence of a saddle point at o = 1/2V suggests 
the most convenient equivalent contour on which to evaluate 
>: If Co is deformed so that it coincides with the 


steepest descents path from the saddle point, then divergent 


it 





integrals are avoided as the integrand of I, becomes 
exponentially small on such a path (granting suitable 


mecerictions on £(o)). This Contour domshoewne> le. = 


“02 o-plane 
Vo <n 
I 
: C 
: R 
I 
1 
nn nl gr 0 
1 
-j 
‘ 
o=1/2V 
| 
I 
ae “5 3 
10, o-plane 
3 ee 2 
| 
| 
J 
{ 
{ 
I 
1 
all oo 
al 








Portions of Co that join the steepest descents path to 
the origin and to the positive real axis at infinity have 


been labeled Ci and Cp: The location of the saddle point 


relative to the simple pole at o = 1/2V depends on the 


parameter V = x/t. Therefore, two general cases are shown. 
The closing contours Ci and Ce are chosen in the fol- 


lowing way: 


il) Cy is the steepest descents path from the origin to 


infinity and intersects the contour Cy at Tei derreteny 


2) Ce isthe arc Of #:he Crecle Ron jOining the contour 


Cc, to the positive real axis. 


Consider the case V = x/t < 1/2. Application of 
Cauchy's integral theorem and the residue theorem lead to 


the following expression: 


-~ itRes(l1) = Lo + Le + Le - 271i1Res(l1) - 271Res (~i) 
1 Z R 


Roo 


(9) I, 


The asymptotic properties of each term in the above 
equality can now be evaluated. 
=t ~-itv 
e e 


Res(-1) = 


lim a Pi Toyo) mle 


+1 
_o* =-1 


This term is exponentially small and can be neglected 


in comparison to terms algebraically small. 





nes(1) = tim, [ (o-])o? jit(o?v-o) ] _ et UF) 
ot) eee —— 


f 


> Q as R>o 
CR 


It is easily shown (see Appendix C) that J. = OiCE a5) 
J 
and is therefore small compared with other contributions to 


I.- Tnerefore, (9) reduces to 


(10) I. - | = An ei (x-t) + 0(t7?) 
&y 


Consider 


(11) | One eit(o*V-o) 4. 
C5 ot=-1 


where V = x/t is a parameter to be specified and is less 
than 1/2 for the case under consideration. Introducing the 
change in variable - i + it = o’°V -o, the following con- 


Siderations apply: 


ie. recs 
qy 7 it = ov O 
] tl 
p= me [le Vivi? } 5a, 124 yt? 1, 7 = V4 
8 
at AV 





On C., in the upper half plane and in the half 


2 
plane to the right of C5 


l do Y 
oS eee a 
Avian 








On Cc. in the lower half plane and in the half plane 
to the left of C. 
i do _ =, 
g=s7l1l-yt?]=o_, a es 
ag AVG 7s 


Therefore, the integral along C. 1s given by 


2 Bele) oe Oe ete (Clete) 


(12) | — ao = | an 
Co o'’-l C. o’-l 








= Z 2 

wo 6 do mo do 

=e av | a ey | ae tae 
O say at 00 


If the Laurent expansion about T = 0 is obtained, the 
above integrals may be evaluated by integrating term by 
term. It 1s more convenient, however, to subtract off the 
pole at o = 1 and to then expand the remainder of the ex- 
pression about t = 0 so that the effect of the pole may be 
considered separately with a view towards obtaining a uni- 
formly valid result. These expansions are developed in 


Appendix D and are given below. 





! cae -y/8 1 aa ee Y -Y, 
+ + IE T 
—— = —— | tT” + + — + 
g 











Y 2 


el 
(-) {1-4vV-4V7}1t /2 r y{ (1+4V7) *7-8V (1+2VvV) 7} z 
(1+2V) (14+4V7) (1+2V) 27 (1+4v2) 2 


(=) v2 (144V2) 2-4v (142) 2 (3-4V2) }t 2 


im + .... 
(1+2V)? (14+4v?2) ? 


14 = 
its o* do_  yY/g ly =o viaiae Y way); 
— = —— | t* - —— +—-]| — + 
of-1 dt 1-2Vv v 8 1-2V 








fr (l=ay—-4y—e) 72 : y{ (1+4v2) 2-8v (1+2v) 7} i 
(1+2V) (1+4V2) (1+2V) 7 (1+4V7) ? 


ee ea Meleted Vid een VAM Vicg) 1 
(1+2v)* (1+4v") = 


Noting that terms containing integral powers of T 
eancel identicad@y Upon substitutienm into (12), the 


resulting expression becomes 








7° oO  =tT 
oc Ltte vee) Bey / See ay eae) da 
ee a 1-2V _f0 1 /*p- = 
ry 
és é a i 
FE = eT at q 4V 2 y¥/8 T 72 v¥/8[1-4V-4V?] 
ha a Soe a? os art |. 
Orne 2 - em O 1-2V (1+2V) (1+4V ) 


3 3 =i 3 ~ 
_ 2/81 (1+4v") 2-4v(2+2V) 2 (3-4V7) 1 A o(1/2) | aoa 
(1+2v) ?(1+4v7) ? 


The first bracketed term in (15) includes the effect 
of the pole at o = 1. The second bracketed term can be 
expressed as a series of Gamma functions and is not in- 
fluenced by the pole at o = 1. Introducing the change of 
72 


variable Tt = s, the integrals in the first bracketed term 


above are transformed to the following: 








_ r ear | ° e Tar _ 
a - = 
: Tote 1-2V a ele 
Y 
— 2 CO — 2 
_ e tSosds am e = 2sds 
O L=2V L=-2V 
s + O = 
x af 
Let s = -s in the second integral of (16). 


Then 





os + 1-2V 





Of, 24 2 ee a. Gowen Ole 
| e “5 2sds i | eo 2sds _ e “5 2sds i | e “© 2sds 
O 




















With the change in variable y = ts the first inte- 
gral in (17) is readily evaluated and the second placed 


in a more convenient form. Equation (17) then reduces to 





1 
(18) 2 pve ee BEY wiz) | 
a Y 
aa e a 
where W(z) = T | aa [I (2) > 0] 
i 
7 (ir2V) 
Y 


W(z) can be represented by the following expression’: 
-z* 
W(z) =e erfc(~iz) 


where erfc(z) = 


eiereamowitz, M. and Stegun, I. A. ed., Handbook of Math. 


freer ons, New York, (1965), p. 297. 





On substituting for W(z), (18) can be placed in the 





form 
1 
(20) 2 | va ee ae ne) j a 
cy 
= 1 ws ai ats 1 
=D | Jj ae in (2%) ¢ ¥ ae (Pe) | 


Therefore, (11) can be expressed as follows for the 


case V = x/t < 1/2: 





fez 1) 2 ; 247 
| ; aitlo V O) a, A 
Co o'<-1 
-~it = (aes) 24 1 ae 
ws e | avait” 2 ani (=) e : erto(i (2% e/*} | + 


: -7/> 
-it/4v [ y/4 1 -Y /4T (1/2) {1-4v-4v7}t 
ie ie ee 


_3 
, rida (1+4v?) 3-4v (1+2v) 2 (3-4v?) HP (3/2) ¢ /2 ef Mee ; 


(27) etd =) 


Substitution of (21) into (10) leads to the following 


expression for (6): 





s 1=2V) 24 
max, t) © Int eS) le ae i erfeli (42/734 





ae 
Y—{1- 4v-4v?} 
erraoy ti eiees P(1/2)t- 


tee T 


oa 7? Pa 


3 VAS ee 3 es 
+ 2iv 2L(it4av*) *-4v t2v) * (3-4v") I}, (3 3) £7 Fon 01 ey |} (veh) 
(1+2v) > (1+4v?) ° 


where y = VY4iV = Y4 x/t 


Consider the case V = x/t > 1/2. Application of 


Cauchy's integral theorem and the residue theorem leads to 


+] + | 
<5 CR 


R->0o 


(23) I. - int Res(l) = | 


a 


Similar considerations regarding | and | 
a CR 
apply as for the case V = x/t < 1/2. Therefore, the 


asymptotic behavior of I, for large t or x is given by 


Zn 


2 2 asi i 
(24) = | aCe Sema OR”) +0 (om) 
Cane =1 : 

2 


Proceeding as for the case V < 1/2, introduce the 
a; te. aS 
change in variable - q +t it =o VO ment prevLOuUS (COs 


clusions regarding the branches and branch cuts ofo(t) 


apply, hence 





(25) 





2 ; 2 
| O aitlo V O) ay _ 
C 











== sale 2 2 a 
me oOo elo) co 86 do 

= e 4V | Ve e tT a 4 | — = aE ea 
eee ct o o*=-1 dt 


The Laurent expansions previously obtained about T = 0 


are valid, hence 





(26) 
—-it a) _1 
| ve AV [14s r(1/2)t” 2 _ y/4{i-av-av2}r(1/2) t772 _ 
os 1-2V (1+2V) (1+4V7) 


2 
— y 3/4{(14+4V*) 2-4V (14+2Vv) 2 (3-4V2) } r(3/2)t7 2 rn o(t7 2) _ 
(1+2v) °(1+4V?) ° 


sete 








Sere oc = 2 CO — 2 = 
_ Y/8 S 4V | e "S$ 2sds - | e *S”2sds 
1-2Vv Ree ie ae an 
x ‘ 


The remaining integrals are treated in a similar manner 
as before except that the substitution s = ~-s should be made 
Mathe first imtegral in lieu of the second integral of (26) 


so that the resulting formulation exists. 
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oO CO = Z 
= 2 See ee eS ds 
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ra JT ane (=) wiz) | 


] = 
where z= te (2520) 





(28) ec, oe : 
=e | va mene in(=2%) ¢ Y erfe{-i (hye 72) | 


The evaluation of (25) is now complete and the wave- 


meeront m(x,t) 1s therefore given by 





i 2 — 
(29) 2 ee Een: 


; eee 7 i 
Riot ae I. =- e ax é ‘i erfe(-i (Yt 7} + 


— ieee / 
ee” Pat eee 
(14+2V) (1+4V?) 


i div “{ (144V2) 3-4v(142V) 3 (3-4V2) } P 


(3) 4772 i 
(P27) a dee =) 2 


—>/, 
+ O(t 4) (vee) 


where y = v4iV 


The asymptotic expansions for n(x,t) have been derived 


fem the two Separate cases of V = x/t < 1/2 and V = x/t-l/2. 
That these two expansions are equivalent can be shown by 


exploiting either of the following symmetry relations: 


(30) erfc(-z) = 2 - erfc(z) 


2 
(31) w(-z) = 2e. - w(z) 


+00 ee 
where w(z) = i/t | S /(z2-t)dt (I Za) 


—OO 





Considering equations (ZZ) sane 2.) = fot) (ca, ee 


symmetry relation (30) is most convenient. 


Hiren 50 je. 














= 1 
ye 2.) acto i 2 - ertc{-i (= Y)¢/7} 
SuUbSt LEU tse (s2) santo -(22)" one ob eagmis 
3) : 
ni Xp I, a qi (x-t) 
-it* _,1-2V, 2 
in _ 4x ( y pe /2} 
a1 2 - erfc{- era me we t + 
Bree oy Hue a 
ce) ale cc 2) | (vy< 172) 
Therefore 
= pee 1-2V, 2 
(34) | | : 
mex ee) ie =- e ax e Y erfo{-i (=) a 


aah es 
a Foe) | } (gaa) 


However, (34) is identical to (29) - the case V > 1/2. 
Therefore, the restriction V Z 1/2 is unnecessary and can 


be removed and either expression completely describes the 





waves. The final result then is uniformly valid about 

the wavefront whose position is defined by V = x/t = 1/2. 
Some simplification of n(x t) is possible and there are 
many equivalent expressions for it. A particularly compact 


expression is given below. 


7 


(35) ett 0% Ii = ei (x-t) ertedy <(i-ty ex)e : ] + 


“(ee 7 z) ie t{l-4x/t-4 2742447 72 -*% 
Se a eimai a alia 
2x7 ee 7) 


For the cases V << 1/2, V >> 1/2, and V = 1/2, expres- 
sions for the waveheight n(x,t) can be developed in an al- 
ternative manner. The effect of the simple pole ato = 1 
is either negligible or an easily evaluated residue contri- 
bution in these cases. These alternative expansions for 
n({x,t) have intrinsic merit and in addition provide a check 
on the validity of the uniformly valid expression that has 


been derived. They are available in Appendix E. 





EXPERIMENTAL PROCEDURE 


Experiments were conducted in the Ship Model Towing 
Tank located in Building 48, Hydrodynamics Laboratory. The 
purpose of the experiments was to provide a comparison be- 
tween the waves generated by a wavemaker and those pre- 


dicted by the derived theoretical expression. 


Control of the oscillatory motions of the wavemaker 
paddle was accomplished by utilizing the sinusoid voltage 
output of a function generator as the input to the paddle 
motion control mechanism - a hydraulic device. A voltage 
output proportional to the deflection of the wavemaker 
paddle was fed back and compared to the input voltage pro- 
viding close synchronization of the paddle motion to the 
desired motion. The amplitude of horizontal displacement 
ol the paddle at the water surface was measured. Both 
frequency of paddle oscillations and maximum amplitude of 
displacement were controllable parameters throughout the 
experiment. It was not found necessary to vary the maximum 
amplitude of displacement from its initial value of 1 inch, 
however. Paddle frequency was varied from 0.8 cycles per 
second (cps) to 2.0 cps. Above a frequency of 1.2 cps non- 
linear effects were encountered and severe vibrations of 
the wavemaker machinery noted. Transverse wave patterns 


were also induced indicating that significant three- 





dimensional features were present in the flow. Asa 
result, all experimental data for runs above 1.2 cps was 
not evaluated. Experiments were not conducted for fre- 
quencies below 0.8 cps so as to avoid the effects of 


menite towing tank depth. 


The height of the waves generated by the wavemaker 
was measured with waveprobes. These instruments are var- 
lable resistance devices whose electrical resistance changes 
with depth of immersion into the water. The voltage output 
of the waveprobes was connected to the input of a strip 
chart recorder (Sanborn, series 7700). Thus a continuous 
record of the waves at a particular location from the wave- 
maker paddle versus time is obtained. Calibration of the 
waveprobes was accomplished by immersing the probes a 
measured distance into the water and adjusting the deflec- 
Bwerom of the strip chart recorder scribe to a convenient 
value. The deflection of the wavemaker boundary was also 
converted to a voltage signal and connected to the chart 
recorder input so that the motions of the wavemaker were 


continuously monitored. 





RESULTS 


To facilitate comparison of the waves predicted by 
equation (35) and those measured experimentally, all dimen- 
Sional quantities are normalized to a consistent base and 
the results illustrated graphically. The transformation 
equations relating dimensional to nondimensional quantities 
are those used in the theoretical analysis in the derivation 
of equation (35) and are as follows: 


Ree 
wx 
xX = where w = frequency of oscillation 


J of wavemaker 





acceleration due to gravity 


ae 
I 
= 
| 
Q 
I 


S = amplitude of oscillation of 
wavemaker at the free surface 


= 
S| 





= 

I 
N 
WY 


All bar quantities are dimensional 


Predicted values of the wave height are readily com- 
puted manually using tables to determine values of the error 
Function complement for complex arguments. Equation (35) is 
well suited, however, for computer evaluation of n(x,t). 
Computer programs for this purpose are presented in Appendix 
G. In addition, a function subprogram for the error func- 


tion complement for complex arguments is listed separately. 


Computer facilities were also used to plot the theo- 


retical and experimental values of wave height using the 





33 
Stromberg-Carlson 4020 (SC-4020) cathede ray tubes oreiugc 
device. Graphs of experimental results are computer faired 
with points corresponding to normalized, measured values 
of wave height encircled. Comparison of faired to experi- 
mental results show that the faired results differ very 
little from measured values. Reproductions of the wave 


records obtained during the experiments are available in 


Appendix F. 


Figures 1 through 5 compare theoretical and experimental 
results for five values of x. Wave height n(x,t) is plotted 
versus time t for fixed values of distance x. These figures 
therefore show the timewise distribution of wave height as 
seen by an observer located at a fixed distance from the 
wavemaker. Figure 6 is generated from equation (35) for 
x = 65 and illustrates well all the prominant features of the 
timewise distribution of waves. Figure 7 shows the spacewise 
distribution of wave height for various of time t. These 
figures are arranged sequentially and show the wave front 


propagating out from the wavemaker with increasing time. 
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Figure 6 illustrates well all the prominant features of 


The wave height in- 


the timewise distribution of waves. 


creases monotonically to a maximum value and then enters an 


oscillatory phase until the steady state wave amplitude is 


reached. 


Figure 7 shows the spacewise distribution of wave 


These figures are 


height for various values of time t. 


arranged sequentially and show the wave front propagating 


out into calm water with increasing time. 
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DISCUSSION —— 


The waves predicted by the equation (35) for the wave 
height as a function of distance and time compare well with 
measured values. To an obServer located a fixed distance 
from the wavemaker, the wave height increases monotonically 
with time to a maximum value and then enters an oscillatory 
phase during which the steady state wave amplitude is 
approached. Both measured and predicted wave records show 


this characteristic behavior. 


Measured and predicted values of the wave height cor- 
respond very closely in phase, except for very small values 
of t when equation (35) is no longer valid. However, meas- 
ured values of wave amplitude are consistently larger than 
predicted values. This difference can be attributed to the 
boundary condition imposed on the moving surface of the 
wavemaker in the theoretical development. For convenience 
of solution, this boundary condition was assumed to be the 


mol lowing: 
w*y/g_. 
(1) s=Se Sin(wt) H(t) (=omaoy < 10) 


However, a more realistic approximation to the flow 


near the towing tank wavemaker paddle is given by 


sa BY Sin(wt) H(t) (yen) 
(36) Ss = 
0 (y 


LA 
i 
7 








This is not the exact boundary condition but is a reason- 
able approximation to it. It assumes that displacements of 
the wavemaker from the vertical are small with respect to 
the depth of the wavemaker paddie (approximately equal to 
the depth h of the towing tank wave generator) and is the 


linear approximation to the actual displacement. 


Ee Wavemaker Free Surface 


=—_—_—_ — =— — em a om oh om = — = = sama ———_ = —_——e a —_=—s = oa ee 


h Fluid 


In addition, the theoretical expression for n(x,t) is for 
the case of infinite depth water. The waves generated by 
the towing tank wavemaker have been chosen so that finite 
depth effects nave a negligible effect on the wave features. 
Therefore, the boundary condition given by (36) is a simpli- 


fied model of the physical situation shown above. 


A tneoretical development based on the more realistic 
approximation to the boundary condition poses no additional 
complication to the problem but does significantly add to 


the algebraic tedium and perhaps obscures evaluation of the 





cs 





results somewhat. Of primary concern is the determination 
as to just how the use of boundary condition (1) above in 
lieu of (36) has influenced the problem. Physical consid- 
erations suggest that assuming an exponential decrease in 
Fluid velocity with depth instead of the more realistic 
linear decay should result in wave amplitude predictions of 
smaller magnitude than measured values. As less energy is 
imparted to the fluid near the moving surface (average 
amplitude of oscillation being smaller), wave amplitudes 
must be correspondingly smaller as energy considerations are 
ultimately related to the square of wave amplitude. This is 


substantiated by experimental results. 


Partial verification of the physical argument presented 
above can be readily obtained in spite of the algebraic 
complexities of the more realistic boundary condition. The 
steady state wave amplitude can be easily obtained by let- 
ting t approach infinity and evaluating the remaining terms. 
For the case at hand this amounts to determining the residue 
contributions to the integral representation for the wave 


Merght. 


The solution to the unsteady wavemaker problem (in- 
finite depth water) can be placed in a more general form 
where the functional dependence of the wavemaker paddle 


displacement on depth is left unspecified (subject to the 





restrictions previously mentioned on this function). One 


Gormulatlon Of it is the following: 


(37) 
i (mer , a 
: = on 
, , Ww VG 
, sin(wt)+sin(/gu S| An 
w+ vYgu 
w 
For F(y) =e J ; 
| e "YF (y) So 
O u + w7/g 
and equation (2) is rederived. 
vo y2ch 
For F(y) = 
0 Vaeieie: bh 
| e YF (y) dy = (uh + ges =) 
0 hu2 


Therefore, the wave height n(x,t) is given by 


(38) esp = 
_ 2ws [ funten = 2) cos (ux) | “34 sin(/gu t)-w sin WE) | au 
T 0 hu? gu - w? 





Introducing the change in variable u = g7/g and 
using non-dimensional variables x = w?/g xX; © = Wene = Jae 
fe= w/o h, and n = mn/2s enesobpcaine. 


(39) 00 _o2h 
nix ee Se | {o*hte ~1} ser (O° xtot) 


ho* (o+1) (0-1) 


yee oye a DS 
ee eee ee) (oe) } do 


Following methods similar to those used in deriving 


jo, as t > 
(40) nix) eee La { in Res[1] } 


(h + eh -~ 1) 7 
h 


% cos (x-t) 


where h = depth of towing tank. 


For the measured depth of 43.5 inches (3.62 feet) and 


wavemaker paddle frequency of 0.9787 cps, h = 4.25. 


(ies ee oa 


4.25 


fot) n(x,t) x Tt Ccos(x-t) = 2.4 cos (x-t) 


Equation (41) compares favorably with measured steady 


state values of wave amplitude corresponding to a paddle 





frequency of 0.9787 cps. These experimentally obtained 


values of wave height are given by: 
(42) Nix,t) = 2235- Coc —t) 


The ratio of the steady state wave height associated 
With boundary condition (36) to that associated with 
mBeundary condition (ly is 2.4/1.57 = 1753. Tf this ratio 
is assumed to hold in a very loose sense throughout the 
transient zone, then the predicted values of wave height 


are very close to measured values. 


For large values of t or x, equation (35) for n(x,t) 
must approach the historical solutions of Lamb and of 
Cauchy and Poisson. This property can be shown by expanding 
the erfc for large positive and negative values of its ar- 


gument. Letting t approach ~ one obtains: 
(43) Mx, ©) hem 2aCOs (x—t) 


These are the outgoing waves at infinity that Lamb derived. 


Letting xX approach » one obtains as the first order term: 


sq? -¥ 
(44) Mex ate) ye x “2 cos(gt?/4x - 1/4) 
2wtT 


Equation (44) exhibits the general features of the 


Cauchy/Poisson solution. 





CONCLUSIONS 


In his analysis of the transient gravity wave response 
to an oscillating pressure, Miles? described many of the 
qualitative features of the generated waves as were observed 
in the towing tank and predicted by the derived expression 
for n(x,t), Equation (35). The wave amplitudes associated 
with a concentrated oscillating pressure also show the mono- 
tonic increase to a maximum value followed by an oscillatory 


phase to the steady state amplitude. 


In addition, Miles defines a rise-time T as the time 
for the timewise envelope to rise from 0.1 of its steady 
state value to its first maximum, and a width X of the 
transition zone for the spacewise envelope as the distance 
between the point where the wave height reaches 0.1 of its 
steady state value and its maximum amplitude. The expres- 


sions for T and X are as follows: 


12 


v4 
(45) i 2 eo (x/qne — 


2 


(46) ere ee 


It is of interest to determine if the theory developed for 
the transient wave response to a wavemaker supports these 


observations or modifies them in any meaningful way. 


Miles, Opel. ep pre L450". 


Equation (35) gives the transient wave response to a 


wavemaker and can be rewritten as follows: 


mig (x,t) %» Li (F eta t) ace eae 2a ean + oc /)) 


Terms of order o(t” 72) can be neglected for sufficiently 
large t (or x). It is convenient to treat the remainder of 
(47) as an outgoing plane wave with amplitude and phase 
given functions of x and t. With this view in mind, (47) 


Gan be rewritten as follows: 


(48) n(x,t) * A cos[{(x-t) + yp] 
where 

(49) B/E Bee sale eel a | 

(50) W = Arg { eee i llsie to) eel? 1 } 


Lixpressions corresponding to (45) and (46) can be 
developed for the transient wave response to a wavemaker 
by investigating the behavior of A and using Miles' 


feminictions Of WX and T. 


Thewbehavior Ofek ase a. Eunece on or U =) <(J-G/ 2x) 1s 


shown below. 
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Wave Envelope 


Timewise 


oe 


Piece 


A is equal to 0.1 of its steady state value (A = 0.15708 = 


A reaches its maximum value when 


ttl J Oe 


O0.1.7/2) when U 


the rise-time T and the 


With these values, 


lee Or. 


In 


A 
Max 


=A ASAE) Fee 


(1-t/2x) 


Se 


width of the transient zone X are readily determined. 
terms of x and t, the maximum and 0.1 steady state ampli- 


tudes occur when 
y 


ie) 


for A= U7 2 Ay iss 


oto 


¥vx (1-t/2x) 


(2 ) 





For xX = <<; 
O 
72 
ce = Zee (2) (1.5255) x, 
i, 
te = 25-52) (2. 195 )geee 
es Oo 
Oils. 
Therefore 
iT Se eae Slt /2 
m O.lss . Oo 
1 
(fp 3) a Oe 
Oo 
ROn we = ts 
Oo 


Vx (1-t./2x, ) ee 21 Ota eee 
apes a) ce oma ae ee Ag.iss 
These relations for Xn reduce to 
Se St Get 2) a a ee tO 
ee - {t.+(2.795) “Ix, 1... + ee = 0 
x and Xo 155 ore tnen given by 


{t+ (-1.5255) *} - eo ee Geo 5>) ie 
xX i 


~ 2 


ee re “ Cie ey eee eee eyes 


(esses ; 





Therefore, 


Soe. - x 
Oimics m 


Sao 2.79572t +7. 82 a 1.5255v2t +2.325 


—s 
— 


2 
For sufficiently large to? 


i ] 
(54) X = 3.06 t ae 3° /2 
O @) 


Equations (53) and (54) can be placed in the same form 


as (45) and (46) by introducing dimensional guantities, thus 


obtaining: 

x l 
(55) T = 9(—2)” 

g 

=3 1 
(56) eo Pie 


Ihe theory developed for the transient wave response 
to a wavemaker thus supports the observations of Miles with 
respect to the asymptotic behavior of the wave envelope 
A(x,t). The conclusions regarding the rise-time T and the 
width of the transient zone X apply without qualification 
to the transient waves generated by both oscillating pres- 


Sure sources and moving surfaces. 





RECOMMENDATIONS 


Several subjects for future analysis are recommended 
which are related to the problem considered by this study. 


They include the following: 


l. Develop a uniformly valid expression for the 
transient wave response to a wavemaker similar to equation 
(35, but using the more realistic approximation to the 
flow near the wavemaker paddle given by equation (36) in 
lieu of (1). Compare the waves given by this expression to 


the experimental results obtained in this study. 


2. Consider the problem of the transient gravity wave 
response to an oscillating pressure treated by Miles but 
using the method developed in this study - first obtaining 
an integral representation for the wave height using Green's 


theorem and an appropriate Green's function. 


3. Consider the transient wave response to a wave- 


maker in a fluid with finite depth. 


4. Consider the transient wave response to an oscillat- 


ing pressure ina fluid with finite depth. 





APPENDIX A — 
THE TIME-DEPENDENT GREEN'S FUNCTION 


The solution of the ene problem of the infin- 
itesimal, irrotational, unsteady motion of a fluid with a 
free surface cannot usually be given in explicit form; 
however, it can be shown that the problem can be reduced to 
that of finding an appropriate Green's Function. The 
essential tool used in solving the problem is the time- 
dependent Green's Function, G(xX,y,2;&,n,C0;t,t), which is of 
intrinsic interest itself. It represents the velocity 
potential which is the solution to the following problem: 
At time t = T a disturbance in the nature of a source is 
originated at a point (x,y,z) which coincides with (&,n,7). 
It is assumed that the pressure on the free surface is 
always zero and that no immersed surfaces are present. G 
satisfies homogeneous boundary conditions on all boundaries 
and appropriate conditions at infinity. After G is con- 
structed, the velocity potential for the given problem can 
be represented in terms of initial conditions and boundary 
conditions. If no submerged bodies are present an explicit 
solution is obtained; otherwise the problem is reduced to 
solving an integral equation.° The two-dimensional time- 
dependent Green's Function will now be developed for 


maifinite depth fluid and fluid of finite depth h. 


a ; a 
Finkelstein, A. B., "The Initial Value Problem for 
Transient Water Waves", Comm. on Pure and Applied Math., 


Vol. X%, pope, oli 522. 








G(x,y;5,0,;t7 = Antinite Beet rl iaee 56 


G must satisfy the following governing equations: 


2 2 
G G 00 0 
VG = oF + a SSE) Goi) Po 
ox y nics Se 
(here 6 represents the Dirac Delta Function) 
Gi, + gG = 0 on y = 0, the free surface 
Y boundary condition 
Ga ASS yas 5 = co 
G = os = 0 for t = tT, y = 0, the initial 
conditions 


Introducing the Fourier transform f£f of G 
oo e 
f = | Ge > Ttax 


one obtains 


G= | fer **au 
TY} we 


Then substituting for G its Fourier transform representa- 


tion, the governing equations become 


1) ao - u’f) = 6(y - n) 


2 fee + gf, = 0 on y = 0 
Bye 20 See ea 
4) f =f, = 0 for tt ="t, y = 0 


The general solution to the above problem can be ex- 


pressed in the following alternative ways: 


oe = ae"Y + be WY = a'cosh(uy) + b'sinh (uy) 


a"cosh[u(yth)] + b"sinh[u(yth) ] 





57 
for y < n, solve equation (1) subject) ce equation) only. 


eer = a eltly 


for y > n, solve equation (1) subject to equation (2) only. 


c 
enor = ¢ cosh(uy) - ao Sinn (uy) 

Require f£f to be continuous at y = n and that the derivative 

of f have a discontinuity of magnitude -1 at y = n, that is 
ettSe - eltbs = -] 


From the continuity requirement above 


Cc 


A ae cosh(un) - =e Sinh (un) ] (y < n) 
elté, as 2e 
C 
el4lnf¢ cosh(uy) - a sinh (uy) J (y > n) 


From the requirement on the derivative of f 


re + guc =g 


subject to c=c, = 0 fort. ="{ 


The solution for c is 


CO = Sas emia) 1) 


G [|e 


ec — 
Therefore, the expression for f reduces to 
. Jun 
eluly|- sinh (un) ae ———{1 wt cos [Vga(t-1) 1) | (y<n) 


u 
oe uly 
eftin| some 5 oye, cos [¥ga(t-1)1}| (yon) 


ells ¢ = 





mnverting tO obtain Gix, wee > ca 


rf" etY 
7 + a Sinh (un) cos {[u(x-€) ]du (y < n) 
Cc re gun ts 
- 2) —— Simh iy cos uc] ck (y > n) 
O 
eo eee 
+ +}, vie weOscy Gu(t-T jy) COs six 2) au 


Recalling the known integral 


eee ak) aes = 
| Sn 2° cos (mx ) ex 16g et oe nS 
: sa Ya2 + m2 
then G may be expressed as 
Gay 2 Gee 
CX, Gitte) See ao log vata) = eae eh 
(CHE)? + (yFn)? 


wi [- BN Avig ls eee cos [u(x-EJdu 


Cy oe) eee ue OF Ch ne eee emma. 


G must satisfy the following governing equations: 








2 2 = 
VG = 23 + 2S = 6(x-A6(y-n) for “? ~ ¥ “© 

Ox oe =0 < X < +00 
Gye + gG = 0 on y = 0, the free surface 

¥ boundary condition 
G, = 0 on y = —-h, the bottom boundary 
condi tron 

G= G, = 0 POW = ely) ty = Ole Clie siete au 


conditions 





Substituting for G its Fourier transform representation f, 


the governing equations become 


el ee 
) ( yy ) (y n) 
a fie + oa = 0 on y = 0 
5) ae) on = <-h 
) y y 
4) f= fF, = 0 Lore t= T poy —Ag 


For y < n, the solution to eg. 1) subject to eq. 3) is 
ae 
a f = a cosh{[u(yth) ] 
For y > n, the solution to eg. 1) subject to eg. 2) Is 


c 
ee = c cosh (uy) - = sinh (uy) 


Require f to be continuous at y = n and that the 
Memivative of £ have@avdiscontinuity of magnitude —-1, 


eat is 


Poe = gon ; | 
¥}y=n+ ¥y=n- 


From the continuity requirement above 


Cc 

| cosh{u(yth) }[c cosh(un) - a Sigh (an Ky een 
wut 

e f = 
coshiu(nth) tic cosh(uy) - <= Selena lis, eee Gy a alee 





From the requirement on the derivative of f 


os Oe 
cosh (uh) 


Ch. + gu tanh (uh) c 


Subject to c = Cc, = 0 fe Ga 


il 
-j 


iiae solution for ¢ is 


_ 1 7 is 
Cc = a Sinh(uh) el coslo(t vt alt 
where o = Ygu tanh (uh) 


Therefore, the expression for £ reduces to 


cosh [u (y+h) ] - ee + 
(y < n) 
‘ cosh[u(nth)]{l - cos[o(t-rt)]} 
Us siicmnicosnt ah) 
alts. me 
cosh[u(n+h) ] - Su =t 
> n) 


A COSh Uyec)) tleecOs Nosce= tie 


u sinh (uh) cosh (uh) 


Maverting to obtain G(x,y;2,n;t,T) 


2 [ cosh [u(y+h) ] sinh (un) 
T Jo 


SeeSan an 2OS [0 (x= 6 ocean) 
ol x ts 
2 cosh [u (nth) ]sinh (uy) 
aaa I. oS Sn Sos ut xc) dusty > 1) 


\ 


1 {cosh [u (nth) ]cosh[u(y+h) ] [l-cos{o(t-t) }] a 
+ t| cose tu Quin) Teese ta ly sh) | ~TEenhtany Cos fu lx €)]du 





Recalling the useful identity 
cosh[u(yth)] = e“Ycosh(uh) - a “Aginn (uy) 


then G can be expressed in the following form 


1 {~ e“Ysinh(un) 
5 u ) 


ra os lu(x=e)jdu (y < n) 
a oe i 
= a | SAN eos (a (x—€ jalan (y > n) 
O 
lite e “4. inh (uy) sinh (un) 
+ a a ——7T eoshtuhn) cos [{u(x-&)]Jdu + 
1{ cosh [u(y+h) Jcosh[u (nth) ] . [l-cos{o(t-t) }] : 
+ a Reais Bese Canhtuny oe oo €) jd 


where o = Ygu tanh(uh) 


Using the known integral expression for the log given 


previously, G can finally be expressed as 


(ET gn? 
Gio se eee = - 5 log ae 
(CRETE yen)? 


oS 
ee | e  sinh(un)sinh (uy) Qooru(x-£)]au + 
0 


u cosh (uh) 


i “cosh [u (nth) ]cosh [u ( +h)] [l-cos{o(t-t) }] a 
e 2] Soom ai u tanh (uh) Se 1c 


where o = Ygu tanh (uh) 





APPENDIXGS 


THE UNSTEADY WAVEMAKER SOLUTION 


Consider a semi-infinite body of incompressible, 
inviscid, fluid having a free surface and infinite depth. 
Let the limiting boundary execute infinitesimal displace- 
ments from an initial vertical plane position given by 
s = SF(y,t) where S has dimensions of length and therefore 
F(y,t) is dimensionless. Initially F(y,0+) = 0. It is 
assumed that F(y,t) is sufficiently smooth and vanishes 
sufficiently rapidly as y approaches -». The problem is 
described mathematically in terms of a velocity potential 


as follows: 


subject to $6 + gd. = 0 ony = 0, the free surface 
Ete Y a 
boundary condition 
a) > O asus ao 
eS SE Gye) on x = 0 
ey = >, = 0 for t = O+; y = 0, the initial 


conditions 
The relevant time-dependent Green's Function asso- 
ciated with the above problem is easily obtained by the 


method of images and is given by: 


Msc ) aaa 7 eee 
Stee 775,07 tpT) =e oo log Ee E) piven (4:5 sae vie) I 


¥(x-&) 2+ (yt+n) 77 (xt) 24 (ytn) ? 


[cos{u(x-&) } + 


m a au (ytn) [1 - cos{Vgu(t-t) }] 
lie u 


+ cos{u(xt+&)}J)du 


Note: The only singularity in the region occupied by the 
Pbupe os ear (ky) Ce jee) lother singular 
points are outside the boundary enclosing the fluid. 

Applying Green's Theorem of vector analysis (in 
2-dimensional form) to the region occupied by the fluid 


exclusive of the singular point (&§,n) one obtains: 


Green's Theorem 


I (6,976. - $,V7,)daA = 
R 


)ds 


35 34 
p (38 ear 





met d4 


G(x eale to) 
Then substituting into the above expression for Green's 


Theorem, 


I] pCV' Opa - ° Oe eo Q (Go, - 6,G,) ds 
im C ae 


where R' = R less the singular point (&,n); Ce is a circle of 


vanishingly small radius enclosing (&,n). 





However, V*¢ = 0 implies V7o, = 0 and ¥-G = Ovexcespreat 


fen) “Thererore 


6 Ser . >, G,)ds a 6 (Gor a $,.G,)ds 
S Cc 
€ 


owe 


On Coy ds = rd@, d/on = - ae ete r = V(x-&)2 + (y =n)? 


With this substitution the above expression reduces to 


: 2T 
TH 
<n —- >,G,) ac iG | | ~(G>,, - oG.) rdég | 


dy 1s regular throughout R; therefore, the only contribu- 
tions to the above limit must come from G and Gy (assuming 


the limit of the integral is equal to the integral of the 


limit) and are given by 


(1) me Gina [| x [=4g2-=] | a0 


ra0 


: [Bei er eee a ; eee 
eee eng . | ae Zi 
Therefore 
} (Go, - $,G,)ds = ———— [ dé = -$,(E,nit) 


Cc 


Thus 


>, (E,n7t) = 6 (¢, 6, - Go, ,) ds 





Integrating the above expression from 0 to T 
T ic 
| obe (eenstvat ~ I, Lo (eS, 7. Sop) 4s | dt 


Note that from the behavior of G and G. at ~ and from the 


prescribed boundary conditions it is easily shown that é 
c 


reduces to 


O O 
ei i 
Cc y=-@2 X=+00 y=0 <= 


On x=0, 3/dn = -d/0x. On y=0, d/dn = +d/dy. Therefore 


18 
| >, (E,nFT) = Gey ee) eri on) = 
O 


— OO 


T ao 
I, I} (eS, 7 Chey oe : I, 6 ler - Coe ay | ene 


Il 
ll 














(y= 0) (x=0) 
oo Let = = 
5 I. | (rey - Daal oe |, | (Shey - ree aye 
(y=0) (x=0 } 
. LS, . : ae 
~ a —— + 
| I Sos i yale [ a t+ 5 ) dt | dx 
(y=0) '- 
ts T zs 
“- : te, eer es | ee, - fo, 8 
(x=0) 





For the problem at hand 


i] 


GL (O,yiE nrt,T) Gy (O,yrésnit,t) = 0 
o, (9,y70) = 0 

th + E> = 0 ony = 0 

G{xX;0; estat) = G,(x,0;7E,nit,T) = 0 


Oi xaC G0 }ea— >, (x,070) = Q 


Y 
o(x,y;0) = 0 (as the fluid is assumed initially 
at rest) 
Therefore 
a8 a 
O(E,n7T) = | [G(O,yrE net, tT) >, (O,yit) - | G, o dt] dy 
0 O 
(x=0) 
Let 
s = § eo” sin(wt) H(t) 
a UoyY = 3 
V = WS emoe COs(wtl te) = o, (O,yit) 
a = -w’s e“oY sin(wt) H(t) + 
+ wS eoY cos (wt)é(t) = bry (Oryit) 
where 
ae Way 


H(t) is the Heaviside step function 


DE ec tiem Dr Tac= De lta funce1on 





Substituting these expressions into ¢(&,n;1T) above, it is 


easily shown that 


-_, _ “wS cos(wt)H(t) (7 _u.y See 
eS eee ean ” e 0" log aE dy + 


S co Un = 
, r2ws | e cos (us) ; | oq (utu.)y ay | . 
Ih 0 u O 


= = 
x | | Ygu sin[Ygu(t-t)] cos(wt)H(t)dt | du 
O = 


Of more direct interest, however, are the waves 
in lieu of the velocity potential. The expression for 
the waveheight is obtained from the linearized free surface 


Thus 


00 * 
mie, t) = ens. | — | | cos{VGu(t-t) }cos (wt) H(t) at | du 
O Oo 0 eo 


or 


_ wsf{ cos (us) [sin(wt)-sin(/gut) 
n(e,= Bf costae) faintwr)-sin gan 


- Sin(wt)+sin(/Ygut) ;du 
o urU, w- ygu w+ /gu 


where u, = w*/g 





APPENDIX se 


ASYMPTOTIC BEHAVIOR OF CERTAIN INTEGRALS 


a. The asymptotic behavior of the integral given below 
1s desired. 


2 ypeeto seem 
| : hla ot) 5. 
Ce Na 


mae contour Cy has been chosen so that it corresponds to 
the steepest descent path from the origin to infinity. 
Introduce the change in variable it = o*(x/t) -o. With 
this substitution, the complex variable T is real along 
this contour. The resulting integral expression is the 


mee LOwlng : 


OO 


| (o?/o*-2) do/dt aaa 


at 

The asymptotic behavior of this integral can be ascertained 
by expanding the integrand about T = 0 and integrating the 
resulting expression term by term. However, the integrand 
is regular about this point and therefore the Taylor series 
expansion is applicable. The series expansion is given 
below and is readily obtained through repeated differentia- 


tion of the integrand. 


O ag Bas 


1 "8 2 
a aia £60) + £°300)7 + Een + 


Oana ee 


The first nonzero term of the Taylor expansion is of 
order 0(1t*). The dominant term in the asymptotic expan- 


sion of | is then proportional to: 
e 
1 





2 se ree a 2 
| . peel edie ae | t2 e “Tar + 
c,° -1 0 


w 


Introducing the change in variable tt = y, the above ex- 


pression can be expressed as follows: 





2 Die: = x 
Ee OO eae) eee y2e Ydy + 

lb 
c,° =a QO 


=—it °T (3) + 


The dominant term in the asymptotic expansion of this 


mmeegral is therefore of order 0(t °). 





Di Consider the following integrals: 
oO . 2 
I, = | £(o) ai lo xtot) do 
O 
_ {f° £0) _ilo?#x+t) 
1, = I, a do 
_ {" £0) _i(o2x-t) 
I, = [ 5 do 





70 
The asymptotic behavior of these integrals can be obtained 
by deforming the path of integration in the complex o plane 


along a contour with the general features of Cy below. 


o-plane 


on 





A particularly convenient equivalent contour for the evalu- 
ation of Ij I,, and I, consists of the steepest descents 


path from the origin to infinity and then along the arc 


R ei? to the real axis. Such a contour is shown below. 


where Can is the steepest 
descents path 
CL Ret? R7*+ © 





For integrals I, and I, the steepest descent path is 
. LT . 

Simply the ray 0 =r er Te For integral Ty: the steepest 

descent path is not as easily described but no further 


complication is added to the analysis. 


As the integrand is in all cases regular at o = Q, it may 
be expanded ina Taylor series expansion about o = 0O and 
the resulting expression integrated term by term. By in- 
spection, the first two terms of the Taylor series are 


equal to zero for I) - Only the first term is zero for 





integrals I, and I,- The contribution to the integral is 
Bien of the order 0 (t=yeertor Gmeeana I, and of order 


mc *) for integrals I, and I,. 


The asymptotic behavior of integrals Ij, I, and I, can be 


described as below: 


0 
I, ~ in Res(l1) + [+ O(t” 7?) = ee ak FON es ) 
C 
R 
(R>+) 
5 e e 
I, v in Res(1) + [+ ec eee + 0(t72) 
: R 
(R+0 ) 
0 . e ‘ 
I, v in Res(1) + a Ge Ge aia One) 
C 
R 


Che) 





APPENDIX D 


LAURENT EXPANSIONS 


To obtain an asymptotic expansion for | 7p Lew isoecou. 


C 
venient to obtain the Laurent expansion of é 


9" do 
o*-1 dt 


about t = 0 such that the effect of the pole at o = 1 has 
been subtracted off with the remainder of the expression 
regular ato = 1- Algebraic difficulties are minimized if 


o*/o*-1 is first expanded by the method of partial frac- 


tions. It is easily shown that 


By ° _ . pts 
0 ON Eee a i + 4 + a oe + a 
Oa 4 +1 O-1 o+tl o-1 











do, 
- Conswider first, the expansion of ws 70-1] ae 
; = 
O do : 
Ta eal eee a aa A 1 =i + 
Tiel at (l+2iv) + yT a (1-2iv) + er 

+ wy ay ys 
Giee2 ee ye (l-2v) + yt%? Sa 


This expression is easily placed in the form 

















2 
Og do Es l _ if rs 
= + » y/8 | ee eV) | + 

gt-1 dt 1-av ’ 

1 
-h +7 
a ae = + 
8 _ l-2v ees 2 
vf 

















mi % 
= Ea a at 
eee /Ae w/e 1 (ic2v) A meal T /2 E 
Ne |. Y Cee a 
ek 
_ {1-4v - 4v*}t 2 Pepe alien eee), 
(1+2v(1+4v7) (1+2v) 7 (1+4v?) ? 
Fe =y {(1+4v*)2 = Ace eae 72 ae 
(1+2v)* (14+4v?) ? 
do 


Beene expansion Of a/ Ce ES VOD tative cd wine a 9S imuenca, 


manner. 
2 


0 ais», i ae 
— = re Cit + 
moto 


o*-1 dt (l+2iv) - yt (l-2iv) - yt 





= af -Y 

+ + 
), vy, 7) 
2 ell nC gt /2 


(1l+2v) - yT =P) = 


By inspection, the expansion of this function is identical 
to the one previously obtained except for sign changes on 


alternate terms. The resulting expression is readily 





placed in the following form: 74 


2 


ae [ Soe = 
pies 





wee le 








eal at l-2v & Y 
-"/p 2 = "/2 2.2 2 
ad ; = eg ATEN 1s ,» yilltav")* = Bv(itav)*} , 
8 JSON (l4+2v) (1+4v7) (l+2v)? (1+4v?) ? 


i y*{(it4v7)? - Aaa te oe - ne 
(1+2v)* (1+4v7) ? 


Beeeror values of v = x/t << 172 or v = x/t>>1/2 there is no 
particular advantage gained in subtracting off the effect 
of the pole at 0 = 1. As before, algebraic difficulties 
are minimized by expressing 07/o0*-1 as a sum of linear 
terms. 


2 7 _, = 
ot do fi prt yk y 2 Jia 
eo =2 dt oti o-i ot] o-l 


Using the familiar expansion 


l _ = = 
SE et co ee ee 


the following expressions Inay be obtained: 





Z I 1 1 
Sele 3 a /2+ /2 4 
t t+ oye Qit/4 | _t AV Ti tlev") 
seo dt 1-l6v" (1-16vV"*) ? 
+ 


| ee, 
ANAS SRP SSRN EDR eg) 
(1-16v")3 





1 i 1 
a Oo Sea 1 7? " av 2472 (1416v") i 
gi-l dt 1-16v" (1-16v") ? 


: 2\ 2 2\ 4% 9 = 
Qyitstil2 (4-4 (Ave) ote 
(1-16v") ? 


For v = 1/2, the following Laurent expansions apply: 

















2 
O do -2_-1 el 
+ See E T Fee ye le oe 
cieere dt 8 16 32 
3 ln. _¥? uote 3, /2 
cee ice Cee] 1.7 
of do_ a ee -7/ e 
— = - y? Ne ie ct a+ — + 
o*-1 dat 8 16 32 
5 c 37 a 
+ — yl / re u oS = ig ar 
64 L238 256 





APPENDIX E 


ALTERNATIVE EXPANSIONS FOR n (x,t) 


For the cases where V = x/t <<1/2, V >> 1/2, and 
V= 1/2, alternative expressions for the waveheight may be 
developed. The effect of the simple pole at o = 1 is 
either negligible or an easily evaluated residue contribu- 
tion in these cases. These alternative expansions for 
N(x/t) have intrinsic merit and in addition provide a 
check on the uniformly valid expression (35). Equations 
(6), (10) and (24) may be combined and rewritten as follows: 


Met) = 1, {3 ae ee een oe | + 0(¢7>)} (V41/2) 
C 


2 
The asymptotic behavior of n(x,t) can be evaluated for 

V << 1/2 and V >> 1/2 by investigating the behavior of 

Ie, This integral (after a suitable change in variable) 
can be expressed as in (12). Laurent expansions of the 
integrand of (12) subject to the restrictions on V above 
have been developed in section 3 of Appendix D. Substitu- 
tion of these expansions into (12) yields the following 


result: 


2 . 2 3 wes = a) sad /2 
| O aitlo V-9) 3, x ay “26 4V 4 = + 


eyo =i o | (1-16vV") 


4ivV{34+12 (4vV2) *+(4Vv2)"*} 3 : -tT 
Wee aS aE eto (7/4) | er 


Integrating term by term, the following asymptotic 


expansion for jee 1S VObitawmech: 
Z 


(tt _ 7 ee / 
| ~ ay 22 e ay a) aan 
G (1-16Vv") 


aM 
Z 


4iV{3+12(4v2) 2+ (4v2) *} 


=: mes 
r (5) t a4 O(t /2) (VSS1/ 2,V<<172) 
(1-16v*) ? 


tieawavenerght W(x ae mcanemow  Douwhiecen: 


Cx es I an ei (x-t) H(1l-2x/t) + 


a 
, “iG - Pp PrGe ” -Y 
+ (-4)V” e ————— + Q(t “7) WWE) 
(1-16v") 
For the case V = x/t = 1/2, the Laurent expansions 


Given in section 4, Appendix D, apply. The simple pole at 
o = 1 is now coincident with the saddle point at o = 1/2V. 
Therefore, the contour C. must be indented and the residue 


contribution accounted for. The waveheight n(x,t) is 


Given by: 


n(x,t) = Tr aT. qi (xt) _ | Aiea) (v=1/2) 
C 
2 





Substitution of the appropriate Laurent expansions 


imo thei integrand of Wiamey elds: 
Zz 





se T, 7 
Se | 
| <n | fer? + 3v2i fe + he Vo 0 ( (17)] ¢ “tt a. 
Cc 
Z O 
ee i) 3 5 7 
» rr rer (S)t fp ieee eer hee (iene (aye. 12 + o(e” 72) | 
The waveheight n(x,t) can now be expressed as: 
Teens - = at (x-t) i 
-i(e-7) - : 
+ (-)e ; : | 2 cS Oh + O(t. fa) | (Yo v= 17 2) 





APPENDIX F 


EXPERIMENTAL DATA 


During the experiments continuous readings of wave 
height and time were made at two locations in the Ship Model 
Towing Tank with the use of a Sanborn, Series 7700, strip 
Chart recorder.: The motions of the wavemaker paddle were 
also continuously monitored and recorded simultaneously with 
the wave data. Time was measured from the instant of 


initial motion of the wavemaker paddle. 


Reproductions of the wave records obtained as the 
output of the strip chart recorder are presented in order of 
increasing wavemaker frequency. ‘Time 1s measured horizon- 
tally and wave height vertically according to the scales 


indicated. 








—- pe Set eke eer Me ee + 
wee —_ — . —. ae o — apeaannr b= 
ee ree 
} i 
shee eee 
- tae ahr eg Sp | 
Lee at = a a ee iL pes Saar | = ‘ 
t + = bons 
ipo =, al code ee = mes. ee, 





' 
et a4 4 S } - oh apt } Ss <——-— psp fl —! aod pry 
por OY aaa ieee |e ona 
= ee | | eueeeen ee j 
p> ee L ete |i ae = = ri 
ry - ss ~ 1 ; T a ae $a tt 
aces, | co 
~ re tj — eee ae | ra | 
aa ree Tee 
ale es eee 
} 
<> eee Se cI 
jae POSS posse 
2 ees —- a 
= Ss So ee 
‘Leos... Se ee 
a | 3 ee eee 
a 3 SC See 





ais, 2.8 shee pect} es eS 
ee Jeo eee 
<a Pe esa | 
ee a a eer eer 
Spl eee ee a 
, 7, # Se ee ee 
ae Eimer] |) | arate etd 

Sar ei 

eet 

mene 
Sar 
Er 


{ 
4 


AT, 
| 


f 
| 
| 





ABTS) caasucicontrms comeusetion 
' 


BREcoag: 


—— rr 


| 
| 
veh 
r 
| 





ia 
He 


t ane 
+ = ofeee SF eS 
ey § let 
1--4 | ass) ae ated 7 4 

im. 
fee cl + 4 
| Sa ees od 
ot 
e es - j » 
> vo of | 
ws ¢ 
ye gS n | 
~A Nw We scl 5 i 
wo 0 eft ; ' 
0 we 
Oe a 4 } 
e et ' 1% 
r4 au) 
4 a 1 g a 
a a Qo «q 
is oO Fw ‘e — 
~YN “A ae N 
2 t(- '- » 5 * 
> ® eo oO . * 
tn Net q by ‘Ss 
oa "eo me 53 5 
” a 4 
4 ee a a f4 Qe 
2 ° ° a e 
oC ba re “A & 
Oe Oy f+ = 
S < 
> z 
$ 





BOT PA rH 


$ 


weil 





ne - 
. 
I 


1 tt: 
pa 





eles 
fo |e 





| z 
we} 
* 


OD S10eaNOD{SHowsD [aa yas 6 yessN B 








°908 [T = WH OT tetwos saty 


(°Oy [T = wm OT reTeos) 
°43 2O9T°SE = X = 28 aqozy 


(0 T = aw ¢ toteos) 
“33 299T°SE = X = T# aqozg 





ae ee 7 
See: te f a a s 
+ So pas et =" 
cs a 2 
seer eee e 
, ~-tr it) ge sqora 
t- +— -4 rade He 
; 1 edo £926 °o - 
{ if ms em ls = : ari 
| We. oe 3 =§_ ¥ = 
A NY A Vi. alte) 
PIVIViVi VY Py vr a 
H ee fee sorzoa uN 
qr Ee. 
Clee 
t | t i 
wr {ot omer) + yay 
4 ee ie 
Ty sory 7, 
ae ae ae peer d ay i 5 









ComTae mk 


ve 


ei 


11.6667 ft. 
1 ca) 


(sonle: 5 sa = 


Probe #l - X = 


Prebe #2 - X= 38.1667 ft. 


(scale: 10 mm = 1 in.) 












_ fw 
~» 
- g@ 
S ; 
® a . 
a ' i 
i: eee 
a i [on 
e { I } 
9 : “emo. otic 
3 ae a 
a =, 
' o 4.._ 
—e ers ce i a ae 
a ee al l 
: rel 
- a oe . . 
a | ee 
3 | { 
Ee L tant 
a — Ps 
> ‘ ‘ 
z e 
z peep 8 
So : 
& & 
= : 
= = a 
= —— 
Pa } 





PT he rt 
oa Soe Re Leet | 
ee ea 
ee er 

Soe tit} tT Ppa Se 

: fabaperee | ee 
asa anes Be IS Pee eT 
eee eer a 1 eae 
. a |_| ee To ee | | 
— Ha : es ea lia ca 
Hs Ez aS ere 
anaen— mae Tt) Se 

| aa Peer 

HF Bee re 
a Ce eee 

mee | la eT 

ce 


j 
rey ss Sa | 





a 








BUFFS O tii ve YOO 
y 

| 

| 

j 

jim 
ft 
Ji 
| 
J 
4 4 
Vel 
1 | 
| 
ij 
L 
ae 
ike 








cad or a 
BSe% at ss al ca ; ce i 2 | Us oe ae - =| | 
eeetpeeee sated eee egy ana ee tenet } — ie) eee | = — eid eee ee ee eo 
Spt po = __ L 
z plot | oar - a S qt 
& | - : 1} me Ss 3 e bas = 
; Lae > Ee > 
CTT ee = ees | ees 
is a a ee 
2 ai oS RSS eee 2 Te See 
1a! eet >) | eee es 
ie et no ae eae || 
tla [ei re | 
-—.aene ee EEL ft 3 
jee edo eee ae 
Ct ee a oe leer | 
BP pote =coneemer scenes Seen. ase 
Sl eee eo pre 
Pao | am 
pias! 
hie 
rate tt 
Pitta gt 
cp pe 
ZEES 5 
BS se oe 
Pee! tl 
ee LO Z 
ehatet | toe} | ae 
Epa a | sides aa se 7 mrp id — 
eee ft i a eg 
Pe CC eee eee pe 





fee eee ee ger ae 


M 


i 
U 


f' 


iY 
wWaWawa ri 


V 


. 
iS 





a 
a 
wee 
PCS. fs 
err | 
aie 






h. Ly 
MyM, 





e e 
TT eS 
— $8 $7 
s+ Sees 
: ae e ae | 
‘ { a a Ag 
Pent 
fo ca a Os 
igre man 
; ; : t= 6 = 
a 
i 2 <4 = 4 
7 sf &§ ‘ 
© 
> 
i 


wh POTION 


40. New “cre 























2 es 

== = A 1 _ 
Peele | ee = 
a epee Rept | | 








YA 

\ y 
- 
| 








1, 
i 


i beh) £ 
Peer arty 
PAE TT ff 


AAU 


I 


i 


ar lh ete 


a pa 
PS AA 


Sy 
ev 


IE 
ca 


PU 
iM 
alae 


ca 
Et 
Pl 





SH Ht dg 


4 





Ca 


rae = aes 
ER eeete 
ed wy ferry et i 
cee Ct pee ee | eee 
el (Le Se ee eee 
Poh tole CT eee teeter Se ee ee 
aa PEE | 
me ae CCRC eed 
- HS {2a se 
ear eo CO ESS el eee TE 
aff eT dy la Pee cae 
a eee a= ee CE edt, | 
| at ae Ce Seer EST 
: ime ee lle Pet fattstablel 
7 fa eS a = z Pet] deere. 
: Se) Le Tee Lie. 
2) TE | te er el fe >, aerate tT 
ict ECE oe eee eee ee wy Fase 
a B a a Tt lol Pen 
eee Sa oie es 
raf | t aoe ee pa 
Pee ~ Spiers eee aE 
! 2 
y ; e) - _ - a ae 
i a = pe "om 
| laa Soaeees Ay | 
CoC eS a a ae me ("EL 
(Cae | a eee Ce eae et 
eo Seer senor can 
z | ae Weert eee ta: | || | Tia | et: | 
g eal. | | a ee Peer 
S are eS eet | 
” tee | |) ee (Cl. iia it ae 
Sty its Sie ae i | Sx ee ae 
3 | eae ee 8 ll EE -- 
eI SS) oo aera 
Si eel = ae tt eee ft ee 
re at ost A et > - cguee s teekl i 
; I = ' . I ' } 
E} | $8 $7 a >. Peet 
Bi | ‘© a 3 | Pe - all Pi a ea Mes a liet 
at | Paes See fo: Sees Bo | 
{i ae © 5 {4-4 | ae 5 ae eer 
Yj: on oe heb fit! io ee eee 
= e a ° eae Pt 5 | ee [ . ty ape (es oe 
So rd 
p: Qs: toy Poche 
ig el t= t¢2 @ eet ee afccealits: | yet |. of 
{ > & @ ro : { | | 
2 =% <3 § an a sow = es | ++. cad 
ai a \ a ; eee 3 ! tea 
Se . 2 ee | | 8 eee | 
a | & 2 cect 
OME OE ja  — ae 
@ i: 5 & ae f ' 2 ma ! 
> ae ee: ae Oo : 
| t dee : ’ . — —— 4 Sheet oon 
9 ea 
joao . | 
a 














sain of 
“+ 


te . 
ail 
h ia 


BRETORD AG CHAK ty 
te 
| | -| 


ae 
Tt 


I 

| 

oat 
is 


te be 


— —eseedaw $s 
Loe ee 
eee 
Saeeeeeerg 
' 

— a + = 
¢ o- asaessa 
‘ t 

‘ o-—mt Oe 





~ ‘ u mom, 4 : =: 
aoe eee = te | | a ee se 
i | —_.. | _._ ae ete aa are 
it.) ee i Eek le 
1. 4} Se er 
em - See eet eee 












eee 
eee 


oT acter OT Hh 
= - sae 









o oa 
ay <a ¢ : E 
gl Sa ee ; eet | 
4 — 
Ae ae | 
a Ag ei 
8 85 Sa a= 
~~ ode : 
 e 
.c ve =a 
-- Net ; 
s2 f= me! 
2 ae ae 
ie Bie 





v8 


+G CH-@ 


































Tie at at rq 
PLT ae pot 


am TL ea : <i 
po4- dL} J - ; | 
bast |i | ee i 
Lhd oe eee Pf 
fees pp =f mente oe eae ja : i ’ 
Pies | Ee ae 
as ie ie ae Pe sais fell 
| | | > 
ad . pete Seo 
rd Tad “pee 
3 | le | -4 cae 
heel | Hee Set 
a rs j | Bee c a: ct oR 
a a Hee Saeeeee ase 
bai ; a 4 a ere. St i 
eet | aaa Tate ty St Pe 
[ene | a Ee eh 
te Tt ae {+ EE 
eee, | SST Cee | Cee 
b se ee) a an Ld eee) aa te eee as val lage tak 
oe Lp D ree TT eee See cy 
eee o5r eee |) P| : Eee tal 
esr hs | Hes ee TI Peete eT | 
~ftedalet ey | ta SCC Sg Ce eee ints 
lee. | +a ee aes 
omer + =e eee SS rr 
; | +. t Seslsaet oe 
He ee ee EES sitccoone 
tt. tie t. r faz tt ho PStisrts fis 
aoe oT pee eee 
ictal 
m= | gee SE 558 
= eee ae 
a 
eee 
ee 
ees ee ae 
Se 
Se Sa 
See 
PES de aes 
ene 
Pea 
HE 
roe ae 
(CER eer 
CECT rae ae 
ot Cea 
| a a 
ae | aa < Pp 
CEC CTT eS Cte Co a 
LS CCC eee Cer ee Perea 
at eee a era a 
a eee TT pas Pret 
= i 
a 
fa 





GRECUAOIRE CARTE, cosmecdooomms conPomTON 


iat eS a 

a. Clear _ Ce eee | 

cH Gos ees , DO 

i ala See | | ieee =F 8 alae | 
—¢- na tp = —- g 9g 
a. aS! 1 ogee. 3 11) le 
aoe [aa | eee eo 2 ian seo 

- Se eee Sa EEE ee 

f = [ 

Sl eet | ae = t Be ee eee oe |e a ek eae 
= Ss Cae Pee 
ie *RCRMEl ee tt 
Sie ee ae i ees eee ES ae — - 

Peele oe ei were! ss | _ fees os Joe Ee 

F $2 34a ae eo) eee 
Sie 3 et See = 

te so a 
cot wy Bae! < ena : me 
bun \ ofl Bo S + —t—— * +--+ 
| 7 8 ! ~w st A | 8B tent 
‘tm - 92 8 f a ° i 
ie : : a 4 x i aan 
1 4 @ _~} . t~ ce oa 
; ut ‘ poe a a” Poot - 
oe z a [4 fee! 

| ° 2 rel ae 
; | a § 2 uy ! eet BI 
Bie Pe : | 
{ 


wT 
\ 








APPENDIX G 


COMPUTER PROGRAMS 


The facilities of the Computation Center of the 
Massachusetts Institute of Technology were used extensively 
in the comparison of theoretical to experimental results. 
The following programs were developed to calculate and plot 


the waves predicted by equation (35) for n(x,t): 


1) Timewise distribution of wave height for 


fixed distance. 


2) Spacewise distribution of wave height for 


fixed time. 


3) Function subprogram for the error function 


complement for complex arguments. 
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APPENDIX H 


NUMERICAL ANALYSIS OF n(x,t) 


ie snort Time sSolueron-. 

The asymptotic behavior of equation (5) for the wave- 
height n(x,t) is readily evaluated for large values of t 
(or X). Tnis behavior is described by equation (35). There 
appears to be no such expression for small values of t 
(or x). For this case numerical methods must be used with 
the aid of computational facilities a practical necessity. 
For very small values of t, equation (2) may be evaluated in 
its present form - integrating numerically; e.g., using 
Simpson's rule. As t (or xX) grows in magnitude, however, 
this approach becomes unsatisfactory due to the increasingly 


eeerllatory nature of the integrand of (2). 


feeerogranming Pitficultics and Che Solution. 


The numerical integration difficulties posed by the 
highly oscillatory nature of tne integrand of (2) can be 
avoided by employing the contour integration methods and 
appropriate cnanges of variable used in determing the 
asymptotic pehavior of equation (5). Here tne integrals Ij, 


I,, I, and I, (defined by (5)) are regarded as contour in- 


4 


tegrals in the complex o plane and equivalent contours 


which facilitate their evaluation have been determined. 








Integrals I TS sae I, can therefore be represented 


ee) 8 
by (see Appendix C): 


2 











; 2 
I, = inRes(1) + | ae eitlo V+) a. 
C o*t-1 
SD 
; . 2 
I, = inmRes(1l) + ee | : a = do 
C o*-l 
SD 
-j ‘o? 
I = itmRes({l) t+ e LE | ' a *do 
4 Con o'=-1 
Introducing the change of variable it, = o7Vt+o into I, 
and it, =g- into I, and I,, these integrals can be written: 
00 > a 
I) = i17Res(1l) + | oO" do e 0 at 
0 o'-l qt, : 
I, = imtRes(l) + et 2 do _ e ~*~ ar 
o 
0 o'-l AT 
2 
Ty = itRes(1l) + eit] FIORE Z, a Sat 
0 o*-l at - 
2 


When expressed in these forms, the oscillatory nature of 
the integrand has been removed and in fact it now becomes 


exponentially small as De and a become large. Numerical 





evaluation of these integrals can now be accomplished by 


a variety of methods for small and large t (or x). 


A similar method can be used to evaluate I.. Its 
representation is given by equations (9) and (23) for the 
cases V = x/t < 1/2 and V > 1/2. Considering only equation 


(23), I. can be written: 


2 2 2 aa 
O PaatG V ao + | O paleo O) a. 
go t-1 C, go t-1 


I, = inRes(1) + | 
Sal 


(Yo Sb) 
where Cy is the steepest descent path from the origin to 


infinity and Co includes both steepest descent paths from 


the saddle point to infinity. Introducing the change in 


variable it = o*V-o into Ne eel Claes /aqy) oor pitenm o*V-o into 
\ 1 
[.., I, can be expressed as: 
C5 Zz 
a2 dc ~tT 
I. = imtRes(l) + a eet 


Oo ace, 


= a 2 = 
mer ee | ee ee (y= 1/2) 
o o*-1 dt, 


A similar result can be obtained from equation (9) and the 


case V . 1/2 - the only difference being a change in sign of 





the residue term. Numerical evaluation of I. is now 
readily accomplished. Finally, equation (5) which relates 
integrals Ii: I, I, and I, to the waveheight, is used to 


complete the numerical evaluation of n(x,t). 
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